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We present a theory for electronic and magneto-optical properties of n-type Ini-ajMn^^As magnetic 
alloy semiconductors in a high magnetic field, B \\ z. We use an 8-band Pidgeon-Brown model 
generalized to include the wavevector (k^^) dependence of the electronic states as well as s-d and 
p-d exchange interactions with localized Mn d-electrons. Calculated conduction-band Landau levels 
exhibit effective masses and g factors that are strongly dependent on temperature, magnetic field, 
Mn concentration (x), and fc^. Cyclotron resonance (CR) spectra are computed using Fermi's golden 
rule and compared with ultrahigh-magnetic-field (> 50 T) CR experiments, which show that the 
electron CR peak position is sensitive to x. Detailed comparison between theory and experiment 
allowed us to extract the s-d and p-d exchange parameters, a and /3. We find that not only a 
but also /3 affects the electron mass because of the strong interband coupling in this narrow gap 
semiconductor. In addition, we derive analytical expressions for effective masses and g factors 
within the 8-band model. Results indicates that (ot — f3) is the crucial parameter that determines 
the exchange interaction correction to the cyclotron masses. These findings should be useful for 
designing novel devices based on ferromagnetic semiconductors. 



PACS numbers: 75.50.Pp, 78.20.Ls, 78.40.Fy 

I. INTRODUCTION 

Recently, there has been much interest in III-V mag- 
netic semiconductors such as InMnAsi and GaMnAs.^ 
The ferromagnetic exchange coupUng between Mn ions 
in these semiconductors is believed to be mediated by 
free holes that are provided by Mn acceptors. They 
become ferromagnetic at low temperatures and high 
enough Mn concentrations. Recent innovative experi- 
ments have demonstrated the feasibiUty of controUing 
ferromagnetism in these systems by tuning the carrier 
density opticalljj^ and electricallyii Understanding their 
electronic, transport, and optical properties is crucial 
for designing novel ferromagnetic semiconductor devices 
with high Curie temperatures. However, basic band pa- 
rameters such as effective masses and g factors have not 
been accurately determined. InMnAs alloys and their 
heterostructures with AlGaSb, the first grown III-V mag- 
netic semiconductor, serve as a prototype for imple- 
menting electron and hole spin degrees of freedom in 
semiconductor spintronic devices. 

The localized Mn spins strongly influence the dclocal- 
ized conduction and valence band states through the s-d 
and p-d exchange interactions. These interactions are 
usually parameterized as a and (3, respectively»i Deter- 



mining these parameters is important for understanding 
the nature of Mn electron states and their mixing with 
delocalized carrier states. In narrow gap semiconduc- 
tors like InMnAs, due to strong interband mixing, a and 
f3 can influence both the conduction and valence bands. 
This is in contrast to wide gap semiconductors where a 
influences primarily the conduction band and /3 the va- 
lence band. In addition, in narrow gap semiconductors, 
due to the strong interband mixing, the coupling to the 
Mn spins does not affect all Landau levels by the same 
amount. As a result, the electron cyclotron resonance 
(CR) peak can shift as a function of the Mn concentra- 
tion, X. This can be a sensitive method for estimating 
these exchange parameters. For example, a recent CR 
study on Cdi-xMuxT^ showed that the electron mass 
is strongly affected by sp-d hybridization. 

In a recent Rapid Communication^ we reported the 
first observation of electron CR in n-type Ini_2;Mna;As 
films and described the dependence of cyclotron mass on 
X, ranging from to 12%. We observe that the electron 
CR peak shifts to lower field with increasing x. Midin- 
frared interband absorption spectroscopy revealed no sig- 
nificant x-dependence of the band gap. A detailed com- 
parison of experimental results with calculations based 
on a modified Pidgeon-Brown model allowed us to esti- 
mate a and (3 to be 0.5 eV and —1.0 eV, respectively. 
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In this paper, we describe details of the theoretical 
model and comparison with the experiments. We use 
an 8-band Pidgeon-Brown model, which is generalized 
to include the wavevector (fc^) dependence of the elec- 
tronic states as well as s-d and p-d exchange interactions 
with localized Mn d-electrons. Calculated conduction- 
band Landau levels exhibit effective masses and g fac- 
tors that are strongly dependent on temperature, mag- 
netic field, Mn doping x, and kz- At low temperatures 
and high x, the sign of the g factor is positive and its 
magnitude exceeds 100. CR spectra are computed using 
Fermi's golden rule. We also derive analytical expres- 
sions for effective masses and g factors within the 8-band 
model, which indicates that (a — /3) is the crucial param- 
eter that determines the exchange interaction correction 
to the cyclotron masses. These findings should be useful 
for desiging novel devices based on ferromagnetic semi- 
conductors. 
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II. EXPERIMENT 

We studied four ~2-/im-thick lui^^Mn^As films with 
Manganese concentrations a; = 0, 0.025, 0.050 and 0.120 
by ultrahigh-ficld magneto-absorption spectroscopy. The 
films were grown by low temperature molecular beam 
epitaxy on semi-insulating GaAs substrates at 200° C. 
All the samples were n-type and did not show ferromag- 
netism down to 1.5 K. The electron densities and mobil- 
ities deduced from Hall measurements are listed in Table 
in together with the electron cyclotron masses obtained 
at a photon energy of 117 meV (or a wavelength of 10.6 
/Ltm). 

TABLE I: Densities, mobilities, and cyclotron masses for the 
four samples studied. The densities and mobilities are in units 
of cm~^ and cm^/Vs, respectively. The masses were obtained 
at a photon energy of 117 meV (or A — 10.6 ^m). 



Mn content x 





0.025 


0.050 


0.120 


Density (4.2 K) 


~1.0xlO^^ 


1.0x10^** 


0.9x10" 


l.OxlO" 


Density (290 K) 


-1.0x10^^ 


2.1x10^'^ 


1.8x10^^ 


7.0x10" 


Mobility (4.2 K) 


~4000 


1300 


1200 


450 


Mobility (290 K) 


~4000 


400 


375 


450 


m/mo (30 K) 


0.0342 


0.0303 


0.0274 


0.0263 


m/mo (290 K) 


0.0341 


0.0334 


0.0325 


0.0272 



The single-turn coil techniquc^° was used to generate 
ultrahigh magnetic fields with a pulse duration of '^l 
fis. The sample and pick-up coil were placed in a liq- 
uid helium flow cryostat. The single-turn coil breaks in 
the outward direction, leaving the sample, the pick-up 
coil, and the cryostat intact, which allows us to repeat 
such destructive pulsed measurements on the same sam- 
ple. We used the 10.6 fim line from a CO2 laser and 
produced circular polarization using a CdS quarter-wave 
plate. The transmitted radiation was detected by a fast 



FIG. 1: Experimental CR spectra for different Mn concentra- 
tions, X, taken at 290 (a) and 30 K (b). The wavelength of 
the laser was fixed at 10.6 /xm with electron-active circular 
polarization while the magnetic field B was swept. The res- 
onance position shifts to lower B with increasing x. The x 
values are 0%, 2.5%, 5% and 12%. 



HgCdTe photovoltaic detector. Signals from the detector 
and pick-up coil were transmitted via optical fiber to a 
multi-channel digitizer located in a shielded measurement 
room. 

Typical measured CR spectra at 30 K and 290 K are 
shown in Figs. H (a) and (b), respectively. Note that to 
compare the transmission with absorption calculations, 
the transmission increases in the negative y direction. 
Each figure shows spectra for all four samples labelled 
by the corresponding Mn compositions from to 12 %. 
All the samples show pronounced absorption peaks (or 
transmission dips) and the resonance field decreases with 
increasing x. Increasing x from to 12 % results in a 
^--^25 % decrease in cyclotron mass (see Table Pi. It IS 
important to note that at resonance, the densities and 
fields are such that only the lowest Landau level for each 
spin type is occupied (see Fig. |2J|. Thus, all the electrons 
were in the lowest Landau level for a given spin even 
at room temperature, precluding any density-dependent 
mass due to nonparabolicity (expected at zero or low 
magnetic fields) as the cause of the observed trend. 

At high temperatures [e.g.. Fig. ^h)] the x = sample 
clearly shows nonparabolicity-induced CR spin-splitting 
with the weaker (stronger) peak originating from the low- 
est spin-down (spin-up) Landau level, while the other 
three samples do not show such splitting. The reason 
for the absence of splitting in the Mn-doped samples is a 
combination of 1) their low mobilities (which lead to sub- 
stantial broadening) and 2) the large effective g factors 
due to the Mn ions; especially in samples with large x 
only the spin-down level is substantially thermally pop- 
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FIG. 2: Midinfrared transmission spectra for the four samples 
studied, taken at 30 K and 300 K. The shift in band gap for 
the X = sample can be attributed to a pronounced Burstein- 
Moss shift resulting from the large electron density in the 
a; = sample. 



ulated (cf. Fig. 13. 

We also performed midinfrared interband absorption 
measurements at various temperatures using Fourier- 
transform infrared spectroscopy to determine how the 
band gap changes with Mn doping x. In Fig. |21we show 
transmission spectra for the four samples at two tem- 
peratures (30 K and 300 K). While a shift can be seen 
in the bandgap in the x = samples, we attribute this 
to a pronounced Burstein-Moss shift resulting from the 
large electron density of the a; = sample. (As one in- 
crease Mn doping, the electron density diminishes rapidly 
since the Mn ions act as acceptors). The Fermi energy 
for an electron density of 1 x 10^^ cm^"^ with a mass of 
0.023mo is 135 meV (without including the nonparabol- 
icity). This roughly accounts for the observed large shift 
in the absorption edge. Thus, our experimental data in- 
dicates that band gap does not depend strongly on the 
Manganese concentration, x. 



III. THEORY 

In this section we describe our effective mass the- 
ory of the electronic and optical properties of n-type 
Ini_j;Mnj;As alloys. Our method is based on the 
Pidgeon-Brownii effective mass model of narrow gap 
semiconductors in a magnetic field which includes the 
conduction electrons, heavy holes, light holes, and split- 
off holes for a total of 8 bands when spin is taken into 
account. In the original Pidgeon-Brown paper, only the 
zone center kz = states were considered and magnetic 
impurity effects were not included. 

The 8-band Pidgeon-Brown model can be generalized 
to include the wavevector (kz) dependence of the elec- 



tronic states as well as the s-d and p-d exchange interac- 
tions with localized Mn d-electronsi^ 

The magneto-optical and cyclotron absorption can 
then be computed using electronic states obtained from 
the generalized Pidgeon-Brown model and Fermi's golden 
rule. We will describe this in detail below. 



A. Effective mass Hamiltonian 

Following Pigeon and Brown^-'^ we find it convenient to 
separate the 8 Bloch basis states into an upper and lower 
set which decouple at the zone center, i.e. kz — 0- The 
Bloch basis states for the upper set are 
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which correspond to electron spin up, heavy hole spin 
up, light hole spin down, and split off hole spin down. 
Likewise, the Bloch basis states for the lower set are 
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corresponding to electron spin down, heavy hole spin 
down, light hole spin up, and split off hole spin up. 

The effective mass Hamiltonian in bulk zinc blende ma- 
terials is given explicitly in Ref. 13. In the presence of 
a uniform magnetic field B oriented along the z axis, 
the wave vector, k, in the effective mass Hamiltonian is 
replaced by the operator 



-A 



(3) 



where p — ~ih V is the momentum operator and A is 
the vector potential. For the vector potential, we choose 
the Landau gauge 



A = -B y X 



(4) 



from which we obtain B — V x A — B z. 

We introduce the appropriate creation and destruction 
operators 
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71 



(kx + iky) 



(5a) 
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and 



a — ^j-^ {kx iky) . 



The magnetic length, A, is 



(5b) 



S^72 f 2N+l 
mo 2 V A2 



- 2kt 



(12c) 




2mo ^bB 



(6) 




(12d) 



where fj,B = 5.789 x 10~^ eV/Tesla is the Bohr magneton 
and mo is the free electron mass. Using Eqs. (|5a|l and 
(|5bl) to eliminate the operators k^ and ky in the effective 
mass Hamiltonian, we arrive at the Landau Hamiltonian 
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Hr = 



La Lc 



with the submatrices La, i& and Lc given by 
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In Eq. (1121 ). we have neglected a second term in M 
proportional to (72 — 73) (a^)^. We do this for two rea- 
sons: 1) (72 — 73) is small and 2) this term will couple 
different Landau manifolds making it more difficult to 
diagonalize the Hamiltonian. The effect of this term can 
be accounted for later by perturbation theory. 

In Eq. p2|) . the number operator N — a^a. The 
parameters 71, 72, and 73 used here are not the usual 
Luttinger parameters since this is an 8-band model, but 
instead are related to the usual Luttinger parameters 7^ , 



I2 



and through the relations^ 
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(13) 
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(10) 



In Eq. 0, Eg is the bulk band gap, and A is the spin- 
orbit splitting. The Kane momentum matrix element 
V = —j-(S\px\X) is related to the optical matrix param- 
eter, Ep, hyi^ 



73 



QEa 



(15) 



This takes into account the additional coupling of the 
valence bands to the conduction band not present in the 
six band Luttinger model. 

The parameter 74 is related to the conduction band 
electron effective mass, m*, through the relationi^ 
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Er, 



74 



1 



E, 



Ey + A 



(16) 
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Er, 



y mo 2 

The operators A, P, Q, L and M are 



(11) 



Note that 74 can be related to Kane paramater F 

\{S\Px\Ur)\' 



F 



{Ec — Er) 



(17) 



A = 



74 /2A^ + 1 



71 f2N+l 2 



P = — — 



mo 2 V A2 



(12a) 



(12b) 



by 74 = 1 -I- 2F. In Eq. l(T7|) . {S\px\ur) is the momentum 
matrix element between the s-like conduction bands with 
energies near Ec and remote bands r with characteristic 
energies Ej.- The Kane parameter F takes into account 
the higher, remote band contributions to the conduction 
band through second-order perturbation theoryiii 
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The Zeeman Hamiltonian is 



tlz — TTT 



Za 

-Za 



where the 4x4 submatrix Za is given by 



(18) 



The total effective mass Hamiltonian for Ini-a^Mn^jAs 
in a magnetic field directed along the z axis is just the 
sum of the Landau, Zeeman, and sp-d exchange contri- 
butions, i.e. 



H = Hr+Hz + H 



Mn 



(26) 
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(19) 



The value of k used in Eq. (|19|l is related to as defined 
by Luttinger through the relationi^ 



For the Luttinger parameter, 
appr oximat ion^ ^ ■ ^ ^ ' ^ ^ 



73^ 



12' 



7f 



(20) 

we use the 
(21) 



The exchange interaction between the Mn++ d elec- 
trons and the conduction s and valence p electrons is 
treated in the virtual crystal and molecular field approx- 
imation. The resulting Mn exchange Hamiltonian is — 



We note that at — 0, the effective mass Hamiltonian is 
block diagonal with respect to the upper and lower Bloch 
basis sets. 

It is assumed in our calculations that the compensa- 
tion arises from As antisites and hence the effective Mn 
fraction, x, in Eq. H22|l is taken to be equal to the ac- 
tual Mn fraction in the sample. We note that this is 
supported by experimantal evidence showing that InAs 
grown at low temperature (200 °C) is a homogeneous al- 
loy and that the magnetization varies linearly with Mn 
content, xJ^i^ii^ 



B. Energies and wavefunctions 

With the choice of gauge given in eq. , translational 
symmetry in the x direction is broken while translational 
symmetry along the y and z directions is maintained. 
Thus, ky and kz are good quantum numbers and the en- 
velope functions for the effective mass Hamiltonian 12()|) 
can be written as 



Hmti = X No (Sz) 



Da 









(22) 



where x is the Mn concentration, Nq is the number of 
cation sites in the sample, and (Sz) is the average spin 
on a Mn site. The 4x4 submatrix Da is 



Da 



ia 



i/3 



2 






i/3 -^^f3 



(23) 



where a and /3 are the exchange integrals. 

In the paramagnetic phase, the average spin on a Mn 
site is given in the limit of noninteracting spins by 



{Sz) = -SBs{ gS 



^J■BB 

kT 



(24) 



where g = 2 and S' = | for for the 3d^ electrons of the 
Mn++ ion.— The Brillouin function, Bg{x), is defined as 

2S+1 /2S+1 ' 



The antiparallel orientation of B and {Sz) is due to the 
difference in sign of the magnetic moment and the elec- 
tron spin. Since B is directed along the z axis, the aver- 
age Mn spin saturates at {Sz) = — |- 



^i(kyy+k^z) 



A 



an,l,u{kz) 




an,2,u{kz) 




an,3,i^(fcz) 
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an,4,i/(fcz) 


K 
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anfi,u{kz) 




an,7,u{kz) 




an,fi,u{kz) 


i>n-l 



(27) 



In Eq. H27() . n is the Landau quantum number associ- 
ated with the Hamiltonian matrix, 1^ labels the eigenvec- 
tors, A = LxLy is the cross sectional area of the sample 
in the xy plane, <j}n{d are harmonic oscillator eigenfunc- 
tions evaluated at ^ = x — X^ky, and an,m,v(kz) are com- 
plex expansion coefficients for the v-ih. eigenstate which 
depends explicitly on n and kz ■ Note that the wavefunc- 
tions themselves will be given by the envelope functions 
in Eq. |57|with each component multiplied by the corre- 
sponding kz — Q Bloch basis states given in Eqs. and 
0. 

Substituting J-n^i, from Eq. (|27|l into the effective mass 
Schrodinger equation with H given by Eq. H26() . we ob- 
tain a matrix eigenvalue equation. By neglecting the sec- 
ond term in M as described in Eq. H12b ). H is block 
diagonal in the Landau quantum number n. We obtain 
a set of matrix eigenvalue equations 



(28) 
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which can be solved separately for each allowed value 
of the Landau quantum number, n, to obtain the Lan- 
dau levels En^u{kz)- The components of the normalized 
eigenvectors, -Fn,!^, are the expansion coefficients, a^. 

Since the harmonic oscillator functions, 0n'(^), are 
only defined for n' > 0, it follows from Eq. (|27|l that 
Fn^u is defined for n > —1. The energy levels are de- 
noted En^yikz) where n labels the Landau level and v 
labels the eigenenergies belonging to the same Landau 
level in ascending order. 

For n = — 1, we set all coefficients ai to zero except for 
06 in order to prevent harmonic oscillator eigenfunctions 
4>n'iO with n' < from appearing in the wavefunction. 
The eigenfunction in this case is a pure heavy hole spin 
down state and the Hamiltonian is now a 1 x 1 matrix 
whose eigenvalue corresponds to a heavy hole spin down 
Landau level. 

For n = 0, we must set oi = 02 = 07 = as = and the 
Landau levels and envelope functions are then obtained 
by diagonalizing a 4 x 4 Hamiltonian matrix obtained 
by striking out the appropriate rows and columns. For 
71 = 1, the Hamiltonian matrix is 7 x 7 and for n > 2 the 
Hamiltonian matrix is 8 x 8. 

The matrix Hn in Eq. H28|l is the sum of Landau, 
Zeeman, and exchange contributions. The explicit forms 
for the Zeeman and exchange Hamiltonian matrices are 
given in Eqs. (jl8|l and H22|l and are independent of n. 
The explicit form of the Landau Hamiltonian for an ar- 
bitrary value of n is given in Appendix IXI 



C. Analytical solutions using Kane model 

In this section we derive analytical expressions that de- 
scribe the conduction band cyclotron resonance energies 
and ^-factors in Ini_2,Mna;As. 



1. General formalism 

The full Hamiltonian of the problem is given by 
Eq. H2f)|l . The natural way to diagonalize this 8x8 
matrix analytically is to treat the off-diagonal elements 
within perturbation theory. This approach, however, re- 
mains valid only if the off-diagonal elements are much 
smaller than the band gap (we are interested in the con- 
duction band states only). Due to strong s-p coupling 
and the relatively small band gap in InAs this condition 
breaks down quickly with increasing magnetic field. In 
fact, VX~^ exceeds the value of one half of band gap at 
B w 25 T. 

Since we are interested in ultrahigh field cyclotron res- 
onance, we choose another approach similar to that used 
by KaneiSi Namely, we neglect the small terms in the 
Hamiltonian matrix H26|l . arising from the free electron 
kinetic energy and the interaction with remote bands. In 
other words, we neglect all terms proportional to A~^. 
Note that these terms are small compared to the band 



gap. They are proportional to the magnetic field and are 
on the order of - 10 meV at a field of B = 100 T. 

Of course, one of the main drawbacks of the Kane 
Hamiltonian is that the heavy-hole band is flat (neglect- 
ing the free electron terms) and the Luttinger model 
as well as the modified Pidgeon-Brown model discussed 
above, were developed to take into account coupling to 
remote bands thereby giving the heavy hole band the cor- 
rect mass. Since in this work, we are only interested in 
the conduction bands, this approximation is not crucial 
and it will allow us to obtain analytical expressions for 
the conduction band energies. 

Our Hamiltonian then is the sum of two parts. The 
first term is a k • p Hamiltonian that takes into account 
only the interactions between conduction and valence 
bands and the second term is the carrier-magnetic ion 
exchange interaction. Neglecting the second term, solu- 
tions of the Hamiltonian can be found analytically. For 
the second term, even in the limit of saturated (Sz)-, the 
exchange interaction is much smaller than the band gap 
and, thus, it can be treated as a perturbation even in 
high magnetic fields. Therefore, as unperturbed states 
we take the solutions of the Kane-like Hamiltonian and 
consider the Mn s{p)-d exchange interaction to first order 
in perturbation theory. 

Solutions of the Kane Hamiltonian can be written in 
the general form 



E 



1 



fiBBi2n + 1) + ± ^g*iE)iiBB, 



where we have introduced an energy dependent dimen- 
sionles effective mass 



1 



p 



m*{E) 3 
and effective g factor 



1 



E + Ea E + Ea + A 



1 



1 



E + Ea 



E + Ea 



(30) 



(31) 



Note that we have neglected the free electron contri- 
butions to these terms. The upper and lower signs in 
Eq. H29|l correspond to spin up and spin down conduc- 
tion band states, respectively, and the zero of energy is 
chosen to lie at the bottom of the conduction band when 
B = 0. 

Note also, for this section only, we have redefined the 
index n compared to the full model which was discussed 
in the previous section. In Eq. (|29|l for electron spin up 
we follow Pidgeon and Brownii and redefine the Lan- 
dau quantum number n by making the transformation 

— > n -|- 1 so that n — corresponds to the ground 
state Landau level for both spin up and spin down so- 
lutions. While this convention on the numbering of n is 
convenient for discussing the Landau levels in the con- 
duction band, it can be applied only for states in the 
simplified Kane Hamiltonian. For the the more general 
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model where the index n ranges from -1 to oo, the lowest 
spin up conduction Landau level is in the n — 1 manifold 
and the lowest lying spin down level is in the n — man- 
ifold. Thus, our use of the Pidgeon and Brown labelling 



with 



where we introduce the notation x = xNo{Sz) /i. Upper 
and lower signs again correspond to spin up and spin 
down conduction band states respectively. 

Note that from Eq. (|32|l we can immediately gain some 
insight into the effects of the s{p)-d exchange interaction 
on the narrow gap material, i) Since we do not take into 
account variation of the material parameters (such as the 
energy gap) with manganese concentration, x, it follows 
immediately from Eq. (|32|l that we obtain a linear depen- 
dence of the conduction band energies on x. ii) Without 
conduction-valence band mixing, the first order correc- 
tion would be E^^^ = ± xNa as one might expect. (In 
this case, there would be no shift in the cyclotron reso- 
nance energy with Mn doping x since all levels for a given 
spin would be shifted by the same amount). The term 
proportional to /3 is a direct consequence of conduction- 
valence band mixing, iii) It is also seen from Eq. 1)32(1 that 
both band mixing contributions induced by the magnetic 
field, B, and motion in the z-direction (fcz) have the same 
sign and, since a and (3 are of opposite signs, both reduce 

To calculate the conduction band energy spectrum us- 
ing Eqs. (|32|I - H33|) we need to know the energies of the un- 
perturbed problem. In spite of its clear form, Eq. H29|) is 
indeed a third order equation and its general solutions are 



where upper (lower) sign corresponds to electron spin up 
(down) states. In Eq. we have introduced V = 
FA-i = ^HBBEp. 



convention for n is confined to this section only. 

Once solutions of Eq. H29|l are known, the first order 
perturbative corrections due to the carrier-magnetic ion 
s{p)-d exchange interaction is given by 



(32) 



(33) 



quite complicated. Standard approximations usually as- 
sume either strong spin-orbit interaction (A ^ Eg + E)2^ 
or small kinetic energy {E ^ Eg + A).^^ In the present 
paper we are interested in strong magnetic fields where 
E ^ Eg K, Is. and have to use the general solution. How- 
ever, in our situation we can make one more simplifica- 
tion, suitable for the particular case of InAs-based semi- 
conductors. Using the fact that in InAs the energy gap 
is approximately equal to the spin-orbit splitting, we set 
Eg — A. Although this approximation is not necessary to 
obtain analytical solutions, it makes the final expressions 
more readable. 

With these simplifications solutions of Eq. (|29|l . i.e. 
energies of unperturbed states in our model, can be pre- 
sented as 

£;(°)«£;o+T„ (34) 

with 

T. = , , ^ (35) 

The position of the bottom of the Landau subbands have 
the form 



- A, (36) 

I 

The second term in Eq. H34() is the kinetic energy for 
the motion along the magnetic field. It assumes parabolic 
dispersion with the effective masses being different for 



/ E 
= ±xN ia + PtisB-^ 



8n 4 T 5 



(E + Eg)^ 
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N 



En 
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(E + Eg)^ [E + Eg+Af 
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each Landau subband. This approximation is vaHd near 
the bottom of the Landau subbands and is justified when 
<t:Eg + Eo. 

Equations (|32|l - 1)36(1 fully describe, within the 
Kane model, the conduction band energy spectra of 
Ini_^Mn^As in a magnetic field. From the energy spec- 
tra one can easily calculate different physical quantities 
observed in cyclotron resonance experiments, such as cy- 
clotron masses or g factors. However, the general expres- 
sions are rather complicated for direct analysis, but can 
be significantly simplified in the limits of low {V <C A) 
and high {V' ^ A) magnetic fields. First we consider the 
situation at the bottom of the Landau subbands. Then 
we take into account finite fc,. 



2. Low-field and high-field limits 

The cyclotron energies are defined as the splitting be- 
tween the two lowest Landau levels (n = and n = 1) 
and are given by AE = E{n = 1) — E{n = 0). 

The low magnetic field limit of the cyclotron energies 
can be obtained by expanding to second order in powers 
oi V / A 1. The result for spin up states is 



AE 



A-— - - /3)-— - xl3-- 



"3 A2 ^ "^'2 A2 
and for spin down states we obtain 



5 y 2 3 2 

AE^^Al^—+x{a~P), 



9 A2 



3 A2 



2 A2 



x(3 



1 V 



I 2 



9 A2 



(37) 



(38) 



It is seen that in small magnetic fields the cyclotron en- 
ergy splitting increases linearly with the field that corre- 
sponds to simple parabolic dispersions at small fc. 

Likewise, the cyclotron energy in the high magnetic 
field limit is obtained by expanding the energy splitting 
to second order in AjV <C 1. The results are 

Ai^T « y fV3 - 1 + 1 A _ o.24^V^("-/3)0-2l|7' 



for spin up and 



9y' 



V 



V 
(39) 



7+0.03-^)+x(a-/3)0.07A, 



J- ^ 9y 

(40) 

for spin down. At high magnetic fields, nonparabolic- 
ity due to conduction and valence band mixing results 
in a square root dependence of the cyclotron energy on 
magnetic fieldM-S^ 

It is also seen from Eqs. (I37ll - (|40|l that the electron 
energy mainly depends not on two independent exchange 
constants a and /3, but on their difference (a — /3). This 
feature will be discussed in more detail in section HV Dl 

It is also easy to obtain limiting expressions for the spin 
splitting of the conduction band states. For the lowest 



Landau sublevels this spin splitting AEq = E^ {n = Q) — 
Ei{n = 0) is 

A y 2 3^2 1 y 2 

^^°^"y A^+^^""^^""^^2 A^"*^9 A^' ^^^^ 
in the low magnetic field limit and 



AEo^-j(l-^)+x{a + /3)-^x~{a~/3) (42) 



in the high- field limit. The first term in Eq. H41|l cor- 
responds to a well-known low-field negative contribution 
to the electron g factor due to the influence of the spin- 
orbit split valence bands Fg and F7 2'^ and can be ob- 
tained by putting E'g = A in the conventional expression 
of Eq. (|31ll . The second term in Eq. H41|l is opposite in 
sign and describes a contribution due to the s-d exchange 
interaction. This term dominates in low magnetic fields 
and is responsible for the giant spin splitting of the con- 
duction bands observed in diluted magnetic semiconduc- 
tors. With increasing magnetic field (Sz) quickly satu- 
rates and the first term in Eq. (|41|l becomes dominant. 
Under certain conditions this may result in a change of 
the sign of the spin splitting. If we go beyond the low- 
field limit, the first term in Eq. H41|) ceases to be linear 
in the magnetic field and, finally, saturates to —A/3 in 
the high- field limit, as seen in Eq. H42|l . At the same 
time, since exchange constants a and /3 have the oppo- 
site signs, strong valence-conduction band mixing results 
in a reduction of the exchange induced spin splitting of 
the conduction band and a term x{a -\- P) appears in 
Eq. lO instead of 2xa in Eq. |(1T1) . 

Finally, let us note that the spin splitting is different 
for different Landau sublevels. For example, for the first 
excited Landau sublevels (with n = 1) the expressions 
analogous to Eqs. (|41|) and H42|) have the form 



A y 2 9 1/' 2 1 y 2 

AE, « -- — +2xa-xia-f3)- — -x/3-— (43) 



3 A2 

in the low-field limit and 



2 A2 



3 A2 



AEi 



A 
9" 



1 



1 A 



+ xia + (3)+x^^{a~l3) 
.(44) 

in the high- field limit. There are two reasons for this be- 
havior. First, with increasing n the energy separation be- 
tween conduction and valence band states increases and 
this suppresses the influence of the valence band on the 
electron g factor. As a result, in the high-filed limit, the 
corresponding term saturates to —A/9, which is three 
times smaller than for Landau sublevels with n = 0. Sec- 
ond, the Landau quantum number n affects the degree 
of conduction- valence band mixing, resulting in a depen- 
dence of the exchange interaction contribution to the g 
factor on n. 
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3. Finite 

The main physics associated with finite is the pres- 
ence of additional conduction-valence band mixing in- 
duced by free motion along the magnetic field. This ef- 
fect is most pronounced if we consider low- field spin split- 
ting of the conduction band states. Expanding the corre- 
sponding expressions to first order in powers of Tq /A <g; 1 
(Tq is the kinetic energy with conduction band edge 
effective mass) we find that the term 

AE^'^-x^{a-P)^-x^P^ (45) 

should be added to the Eq. (|41() in general case of states 
with kz ^ 0- The last two terms in Eq. I|41|l together 
with Eq. 1451) describe the influence of the s-p coupling 
on the spin splitting of the conduction band Landau sub- 
bands. It is worth noting that both these contributions 
lead to a reduction of the s-d exchange splitting 2xa. Let 
us also note the difference in their manifestation. Mag- 
netic field induced s-p mixing becomes pronounced at rel- 
atively strong magnetic fields. At such fields, however, 
the role of exchange interaction is decreased by the con- 
ventional spin splitting (first term in Eq. 1)41(1 '). At the 
same time, the k^ induced contribution could be impor- 
tant even in small magnetic fields, when the exchange 
interaction plays the dominant role in spin splitting of 
conduction band states. 

In our model, exchange interaction is considered within 
first order perturbation theory. As a result, all physical 
quantities like ij-factors or cyclotron masses are linear 
functions of the manganese concentration x. This means, 
in particular, that if the exchange interaction contribu- 
tion to the spin splitting vanishes for a certain concen- 
tration of magnetic ions, it should vanish at the same kz, 
magnetic field and temperature for any value of x. The 
results of numerical diagonalization of the full Hamilto- 
nian matrix (|26|) . presented in Fig. |S1 supports this con- 
clusion and justifies our model. 



D. Magneto-optical absorption 

In this section, we discuss how we calculate the 
magneto-optical properties. We calculate the magneto- 
optical absorption coefficient at the photon energy huj 
fromS^ 

aihuj) = e2{nu) (46) 

[ncjrir 

where e2{fiuj) is the imaginary part of the dielectric func- 
tion and rir is the index of refraction. The imaginary part 
of the dielectric function is found using Fermi's golden 
rule. The result is 



E 



n.i/;n ,1/ 



dkz |e-P„V (fc.)l^ 



where Ai?";|^^, (fc^) = En'.yikz) - En,u{kz) is the tran- 
sition energy. The function fn,v{kz) in Eq. (|47|) is the 
probability that the state (n, v, kz), with energy £'„,t/(fcz), 
is occupied. It is given by the Fermi distribution function 



l + eM{EnAkz)-Ef)lkTy 



(48) 



The Fermi energy Ef in Eq. ((48|l depends on tem- 
perature and doping. If Nu is the donor concentration 
and Na the acceptor concentration, then the net donor 
concentration, Nq = No — Na, can be either positive 
or negative depending on whether the sample is n- or p- 
type. For a fixed temperature and Fermi level, the net 
donor concentration is 



Nc 



1 



(2^)2A2 ^ J_ 



dkz UnAkz) - (49) 



UnAkz) - fwyikz)) S (Ai?;r,,(fc,) - hu;^ , (47) 



where 5^ ^ 1 if the subband (n, v) is a valence band 
and vanishes if (n, v) is a conduction band. Given the 
net donor concentration and the temperature, the Fermi 
energy can be found from Eq. (|49() using a root finding 
routine. 

Since the envelope functions and vector potential are 
slowly varying over a unit cell, the dominant contribu- 
tions to the optical matrix elements are given by 

^' PnA (^^) ~ E/ ^n,m,fj.ikz) Un' ,rn' ,fj.' (kz) 
m,ra' 

x((/>JV(n,m)|0W(n',m')> (TO|(e-F)|m') (50) 

where e is the unit polarization vector of the radiation, 
o-n,m,i.iikz) are the complex expansion coefficients for the 
envelope functions in Eq. 1(2 7() . and (l)N{n,m) are orthonor- 
malized harmonic oscillator wavefunctions. Their quan- 
tum numbers N(n, m) depend explicitly on n and m as 
defined in Eq. 1(2 7() . In Eq. ((50(1 . we have neglected a 
term that depends on the momentum matrix element, 
('/'jv(n,m)|(e • P)\4'N{n' ,m')) between the oscillator states. 
Owing to strong band mixing in the narrow gap materi- 
als, this term is much smaller than the momentum ma- 
trix elements between the Bloch basis functions, even for 
intraband optical absorption such as for cyclotron reso- 
nance, hence we neglect it. 

The momentum matrix elements {m\Pa:\m'), 
{m\Py\m'), and {m\Pz\m') are the momentum ma- 
trix elements between the Bloch basis functions |m) 
defined in Eqs. and Explicit expressions for 

the momentum matrices Px, Py, and Pz are found in 
Appendix ^ 

In our simulations, we consider e-active circularly po- 
larized light incident along the z axis. For e-active cir- 
cular polarization, the unit polarization vector is e = 
{x — iy) / In performing the integral in ((47() the Dirac 
delta function, 5{x), in Fermi's golden rule is replaced by 
the Lorentzian lineshape function A^(a;) with full width 
at half maximum (FWHM) of 7. 



10 



TABLE II: Material parameters for Ini-^^Mn^^As 



Energy gap (eV) " 

Eg (T = 30 K) 

Eg (T = 290 K) 

Electron effective mass (mo) 

m* 

Luttinger parameters 

it 
it 
it 

Spin-orbit splitting (eV) ' 
A 

Mn s-d and p-d exchange energies (eV) 
A''o a 
No p 

Optical matrix parameter (eV) ' 
Ep 

Refractive index 



0.415 
0.356 

0.022 

20.0 
8.5 
9.2 

7.53 

0.39 

-0.5 
1.0 

21.5 

3.42 



"Eq. ISll with parameters from Ref. |14L 
''Ref.hj 
'^Ref.n 



E. Material parameters 



The material parameters we use are shown in Table Ull 
For the temperature dependent energy gap, Eg, of InAs 
we use the empirical Varshni formulai^iSS 



En — Eq — 



ill 

T + b 



(51) 



with Eo = 0.417 eV, a = 2.76 x lO""* eV/K, and 6 = 93 K. 
As seen in Table I, the experimental cyclotron mass does 
not vary significantly with temperature. Other low field 
CR studies of InAs^S also show a very weak tempera- 
ture dependence of the cyclotron mass. To account for a 
mass that only slightly varies with temperature, we keep 
the mass constant with temperature and adjust the 74 
parameter (Eq. H16|l ) to account for the temperature de- 
pendent band gap. Alternatively, one might wish to keep 
74 constant and vary the k-p optical matrix element Ep. 



IV. RESULTS 

A. Landau levels and g factors 

In Fig. 13 (a), the six lowest lying zone center {k^ = 0) 
Landau conduction subband energies in InAs are shown 
as functions of applied magnetic field. At = 0, the 
submatrix, Lc, in the Landau Hamiltonian, H^, vanishes 
identically and the Hamiltonian, H, in Eq. (|26|l is block 




(Tesia) 



FIG. 3: Zone center Landau conduction subband energies 
at T = 30 K as functions of magnetic field in n-doped 
Ini_^Mn^As for (a) x = 0% and (b) x = 12%. Solid lines 
are spin-up and dashed lines are spin-down levels. The Fermi 
energies are shown as dotted lines for n — lO^'^cm"'' and 



diagonal with respect to the upper and lower Bloch basis 
states in Eqs. (Q) and ||2J). In Fig. Othe spin-up Landau 
levels are shown as solid lines while spin-down Landau 
levels are shown as dashed lines. The Fermi energies for 



tration of n = 10^^ cm"'^, only the lowest Landau sub- 
band is occupied at i? = 60 T, while at n = 10^^ cm"'^, 
the first two Landau subbands are occupied. 

In the absence of Mn impurities, it is well known that 
at low magnetic fields, the spin splittings between the 
electron spin up and spin down Landau levels, at the 
band edge, can be described in terms of an effective gy- 
romagnetic factos^ZiS^ 



9 



= 2 1- 



En 



3Eg Eg 



(52) 



This is just Eq. H31|) with E = Q and the bare g factor 
2 included. The g factor in Eq. (|52|l depends on tem- 
perature through the temperature-dependent band gap. 
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For bulk InAs at T = 30 K we find g* = -14.7 and for 
T — 290 K we have g* = —19.0. At nonzero magnetic 
fields, the expression for the gyromagnetic factor in Eq. 
(|52|l is not correct. As a function of the magnetic field, 
the spin splittings become nonlinear as seen in Fig|3|[a) 
and depend explicitly on the Landau subband level in- 
dices, n. This results from the nonparabolicity in the 
narrow gap material. 

The effect of doping InAs with Mn is shown in Fig. 
Otb) where the Landau subband levels for Ino.88Mno.12As 
are plotted. At low fields, the effect of doping with Mn 
is to alter the gyromagnetic factor. The gyromagnetic 
factor in the presence of Mn impurities ia^ 



9Mn 



X NqU {Sz 



(53) 



Since the sign of a is negative, the gyromagnetic factor 
increases with the applied magnetic field. At sufficiently 
high Mn concentrations and magnetic fields, the low tem- 
perature gyromagnetic factor can become positive and 
the spin splittings reverse as seen in Fig. Ofb). 

We have examined the electron g-factors for the low- 
est spin down and spin up Landau levels with energies 
^0,4(^2) and Ei_Q{kz), respectively. The gyromagnetic 
factor is obtained by diagonalizing the Hamiltonians Ho 
and Hi in the matrix eigenvalue problem H28|) and com- 
puting the gyromagnetic factor as 
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The gyromagnetic factor can depend on the temperature, 
the magnetic field, and the wavevector. 

In Fig. 21 the gyromagnetic factor for the zone cen- 
ter is plotted as a function of the magnetic field for two 
different temperatures and four different values of the 
Mn concentration. In all cases, doping with Mn impu- 
rities is seen to increase the g factor and the effect is 
seen to be sensitive to both the temperature and the Mn 
concentration. The sensitivity to temperature and Mn 
concentration arises from the factor x (Sz) in the Mn 
exchange Hamiltonian (|22|l . In the paramagnetic phase, 
(Sz) is related to the temperature and magnetic field in 
accordance with the simple Brillouin function expression 
ofEq. j21I). 

In Fig. [SI the gyromagnetic factor is plotted as a func- 
tion of the wavevector kz for a magnetic field of B — 10 
T. As in Fig. 01 two temperatures and four Mn concen- 
trations are considered. We note that in Fig.Efb) all the 
curves cross around kz ~ 0.5 nm~^ indicating that the 
exchange interaction contribution to the g factor vanishes 
for this state regardless of the Mn concentration. The 
reason for this curious situation was explained earlier in 
our discussion of the Kane model. 



FIG. 4: Gyromagnetic factors at = for the lowest Lan- 
dau subband in n-type Ini_3;Mna;As as a function of applied 
magnetic field for several values of the Mn concentration, x. 
The upper panel is for T = 30 K and the lower panel is for 
T = 290 K. 



B. Magneto-absorption 

The magneto-absorption spectrum and the bandstruc- 
ture for n-doped Ini_a;Mna;As are shown in the upper and 
lower panels of Fig. for e-active circular polarization. 
In our simulation, the Mn concentration is a; = 12%, the 
external magnetic field B = 34 T, the carrier concentra- 
tion n = 5 X 10^^ cm~^, and the temperature T = 30 K. 
The full width at half maximum (FWHM) is taken to be 4 
meV in our calculation. In the bottom panel of Fig. |2|the 
Landau subbands are plotted as a function of wavevector, 
kz, parallel to the applied magnetic field. The Landau 
subband energies depend only on kz and the bandstruc- 
ture is one dimensional. The Fermi energy is indicated 
by the dotted line and at this carrier concentration and 
temperature, only the lowest lying spin-down conduction 
subband is populated near the zone center. Electrons in 
this partially filled subband can be excited to higher ly- 
ing conduction subbands. A strong An — 1 transition 
(labelled 1) is observed between the filled ground state 
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FIG. 5: Gyromagnetic factors at B = 10 Tesla for the lowest 
Landau subband in n=type Ini-xMn^As as a function of k. 
g factors for several values of Mn concentration are shown. 
The upper panel is for T = 30 K and the lower panel is for 
T = 290 K. 
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FIG. 6: The upper panel shows the absorption spectrum for 
n-doped Ino.88Mno.12As in a magnetic field. The radiation is 
e-circularly polarized and the FWHM is taken to be 4 meV. 
The bandstructure and Fermi energy are shown in the lower 
panel. Three vertical transitions labelled 1, 2, and 3 corre- 
spond to absorption features in the upper panel. Transitions 
between spin up levels are shown as solid lines and transitions 
between spin down levels are shown as dashed lines. 



conduction subband and the first excited spin-down con- 
duction subband. Since these two Landau levels have the 
same curvature, a sharply peaked joint density of states 
results in the sharp peak in magneto-absorption observed 
at a photon energy hui = 0.117 eV. Two valence to con- 
duction absorption features (labelled 2 and 3) are seen 
for hu < 0.5 eV. Transitions between the ground state 
spin-down valence and conduction subbands give rise to 
the feature labelled 2 in the figure. The bandedge for this 
absorption feature depends on the position of the Fermi 
energy due to Fermi blocking effects. Another valence 
to conduction subband transition (labelled 3) is due to 
transitions between the ground state spin-up valence and 
conduction subbands. The valence subband has a char- 
acteristic camel back structure and near the zone center 
has the same curvature as the conduction subband. This 
results in an enhancement in the joint density of states 
near the zone center and gives rise to the peak in the 
absorption spectrum near huj = 0.47 eV. 



C. Cyclotron resonance 

In Fig. 13 we simulate cyclotron resonance experiments 
in n-type InAs for e-active circularly polarized light with 
photon energy ?iuj — 0.117 eV. We assume a temperature 
T = 30 K and a carrier concentration n = 10^^ cm""^. 
The lower panel of Fig. [3shows the four lowest zone cen- 
ter Landau conduction subband energies and the Fermi 
energy as functions of the applied magnetic field. The 
transition at the resonance energy Hlu = 0.117 eV is a 
spin-up An = 1 transition and is indicated by the verti- 
cal line. From the Landau level diagram the resonance 
magnetic field is found to be i? = 34 T. The upper panel 
of Fig. [7|shows the resulting cyclotron resonance absorp- 
tion assuming a FWHM linewidth of 4 meV. There is 
only one resonance line in the cyclotron absorption be- 
cause only the ground state Landau level is occupied at 
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FIG. 7: The upper panel shows the cyclotron absorption as a 
function of magnetic field in n-type InAs. The radiation is e- 
active circularly polarized with huj = 0.117 eV. The FWHM 
linewidth is taken to be 4 meV. The lower panel shows the 
four lowest conduction Landau levels and the Fermi energy as 
a function of applied magnetic field. Solid lines are spin-up 
levels and dashed lines are spin-down. 



low electron densities. 

At higher electron densities multiple hnes can appear 
in the cyclotron absorption as higher lying Landau sub- 
bands are populated. This is the spin-splitting of the 
cyclotron resonance peaks. This is illustrated in Fig. 
|S1 where the electron density, n, has been increased to 
10^* cm~^. We can see in Fig. Othat both the lowest 
lying spin-up and spin-down Landau levels are occupied 
and that two resonance transitions (labelled 1 and 2 in 
Fig. IS)) result. In addition, we also see that the spin- 
down peak which comes in at higher electron densities 
appears to dominate the spin-up peak. This is somewhat 
misleading. The dashed-dotted and dashed lines show 
the contributions to the cylotron resonance absorption 
from transitions 1 and 2 respectively. Both transitions 
have roughly the same strength, but the lineshape for 
transition 1 is highly asymmetric. This results from tak- 
ing into account the contributions to the absorption for 



FIG. 8: Same simulation as in Fig. Qbut with the doping 
density raised to n = 10^**cm~'^. Since the two lowest Landau 
levels are occupied, there are now two cyclotron absorption 
peaks. Peak 1 is the spin up transition and Peak 2 is an ad- 
ditional spin down transition. The dashed-dotted and dashed 
lines show the individual contributions to the cyclotron res- 
onance absorption from transition 1 and 2 respectively. The 
asymmetry of the lineshape results from the fc^ 7^ contribu- 
tions as well as the nonparabolicity. 

kz and also the nonparabolicity of the bands. The 
lineshape for transition 2 is not as assymetric since the 
Fermi energy lies closer to the Landau level edge for this 
transition. Note that if there is substantial scattering so 
that the transition lines are broadended, then the two 
transition will be merged into one and can lead to an 
apparent shift of cyclotron resonance features to higher 
fields. 

The effect on the cyclotron resonance absorption of 
doping InAs with 12% Mn is shown in Fig. |51 Com- 
paring Figs. and IHl we see that the effect of heavily 
doping with Mn while keeping the electron concentra- 
tion, n, fixed is to reverse the spin splitting and hence 
the character and position of the cyclotron absorption 
peak. The cyclotron absorption peak seen in Fig. El is a 
spin-down transition as opposed to a spin-up transition 
and the cyclotron resonance peak is seen to shift from B 



14 




> 

0) 



c 

UJ 



0.7 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 ij-r^ \ \ \ \ \ 

--"-"■"^ ^ _ 
^^^^ ^ 


0.6 














— ^ 

- ^ 




0.5 


— 

^ — " 


Fermi - 


0.4 


~i 1 





20 30 40 50 60 

B (Tesia) 



2 




20 



-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

L(b) A 
= 290 K / 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 U 

X..^ 12% = 




5% = 




V.^ 2.5% = 


— / ' ^ 

— / ' 
^ -L-L-J-atf^l ''l 


0-01 0% i 

1 1 iT-.i-ili-1-i J J i L ^- 



30 40 50 
B (TesIa) 



60 



FIG. 9: Same simulation as in Fig. Qbut with the Manganese 
concentration increased to a; = 12%. The cyclotron absorp- 
tion peak near 32 T is now a spin down An = 1 transition. 



FIG. 10: Cyclotron absorption as a function of magnetic field 
in n-doped Ini_^Mn^As with x = 0%, 2.5%, 5%, and 12% for 
(a) T = 30 K and (b) T = 290 K. The radiation is e-active 
circularly polarized with hui = 0.117 eV. 



= 34 T to B = 32 T. 

In Fig. ^1 we examine more closely the effects of the 
Mn doping concentration, x, on the cyclotron absorption 
spectra of n-doped Ini_a;Mn2,As for e-active circularly 
polarized radiation with photon energy huj — 0.117 eV. 
We plot the cyclotron absorption as a function of applied 
magnetic field for x = 0%, 2.5%, 5%, and 12%. In aU 
cases we assume a fairly narrow FWHM linewidth of 4 
meV. In Fig. ^| (a) we plot the cyclotron absorption 
at T = 30 K. All the solid curves are computed for an 
electron concentration n = 10^^ cm~'^ and are vertically 
offset for clarity. For x = 0%, we plot the cyclotron ab- 
sorption for n ~ 10^^ cm~'^ as a dashed line scaled by 
a factor of 0.01. In Fig. ^| (a), the cyclotron absorp- 
tion curves for x = 0% and x = 12% have already been 
discussed in Figs. El and El As x increases from 0% 
to 12%, the cyclotron absorption peak initially shifts to 
higher fields and is due to a An — 1 spin-up transition 
from the occupied lowest Landau level. At a critical value 
of the Mn concentration, the spin splittings reverse and 
the observed transition is now due to a An = 1 spin- 
down transition. Consequently, the cyclotron absorption 



begins to decrease with increasing Mn concentration. 

To examine the temperature dependence, we also com- 
puted the cyclotron absorption at room temperature for 
several values of x. The curves in Fig. ^1 (b) are the 
same as those in Fig. ^| (a) except that the tempera- 
ture has been increased from T = 30 K to T = 290 K. 
We see that the single cyclotron absorption peaks ob- 
served at low electron density in Fig. ^| (a) split into 
doublets in Fig. llOf b). At room temperature, the two 
lowest Landau levels are thermally populated and both 
the spin-up and spin-down transitions appear in the cy- 
clotron absorption spectra. The shift in the cyclotron 
absorption features with increasing Mn concentration is 
also different at room temperature. As x increases, we 
find that the cyclotron absorption features in Fig. 1101 (h) 
shift to higher fields. As we increase the temperature 
the magnitude of the average Mn spin (Sz) described by 
the Brillouin function in Eq. (|24ll decreases. Since the 
exchange Hamiltonian, Hmu j in Eq. (|22|l is proportional 
to X (Sz), the spin splittings reverse sign at higher Mn 
concentrations. 
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FIG. 11: Calculated electron cyclotron masses for the lowest 
lying spin up and spin down Landau transitions in n-type 
Ini_2:Mni:As at hu) — 0.117 eV as a function of Manganese 
concentration for (a) T = 30 K and (b) T = 290 K. Electron 
cyclotron masses are shown for three sets of a and (3 values. 



D. Electron cyclotron mass 

The electron cyclotron mass, Mcr, for a given cy- 
clotron absorption transition is related to the resonance 
field, B* (cyclotron energy ^bB*), and photon energy, 
hoj, by the definition 

Mcr 2fiBB* 

= — (55) 

Too riLu 

In Fig. 111! the calculated cyclotron masses for the lowest 
spin-down (b set) transitions and spin-up (a set) transi- 
tions are plotted as a function of Mn concentration, x, 
at a photon energy of hui — 0.117 eV. Cyclotron masses 
are computed for several sets of a and /3 values. The 
cyclotron masses in Fig. (a) and (b) correspond to 
the computed cyclotron absorption spectra shown in Fig. 
[TUl (a) and (b) , respectively. In our model, the electron 
cyclotron masses shown depend on the Landau subband 
energies and photon energies and are independent of elec- 
tron concentration (though clearly, the population of a 



given state will depend upon the concentration). 

In the Kane model discussed earlier, the cyclotron en- 
ergy in Ini_2,Mna;As depends linearly on the Mn con- 
centration, X. In addition, the cyclotron mass, which is 
inversely proportional to the cyclotron resonance energy, 
increases for conduction electrons from the (a) subsystem 
and decreases for conduction electrons in the (b) subsys- 
tem. These predictions of the simple Kane model are 
both confirmed as can be seen in Fig. ^2 

The sensitivity of the cyclotron masses to a and /3 
seen in Fig. ^] can be understood if we study the de- 
pendence of the cyclotron energy on the exchange con- 
stants a and (3 in the simple Kane model previously dis- 
cussed. Because InMnAs is a narrow gap semiconductor, 
it is not surprising that the conduction band cyclotron 
energy should depend not only on the conduction band 
exchange constant a but also on the valence band ex- 
change constant /3. This just reflects conduction-valence 
band mixing in the magnetic field. Our Kane model cal- 
culations show, however, that only their difference is im- 
portant and that only one independent constant (a— /3) is 
needed to describe conduction band cyclotron resonance 
with high accuracy. It is seen from Eqs. (|37|l - H4l)|) that 
in the high-field limit the exchange interaction correction 
to the cyclotron energy is proportional to (a — f3) while 
in the low-field limit the term proportional to (a — (3) 
is one order of magnitude larger than the term propor- 
tional to /3. Most of the important corrections to the 
cyclotron energy come from the heavy hole admixture to 
conduction band states. The difference (a — /3) refiects 
the change in the band gap due to the exchange interac- 
tion and affects the degree of this mixing. This suggests 
that the dependence on (a— (3) can be interpretted as just 
a renormalization of the band gap in a magnetic field. 

Note that while Fig. II II shows that the cyclotron peak 
positions are very sensitive to (q — /3) (the figure shows 
how a 20% change in (a — (3) can affect these positions), 
one can only measure the difference a — (3 from these 
measurements and not a and /3 independently. In addi- 
tion, to determine {a — j3) the other parameters (such as 
7i: 72, 73, etc.) must be accurately known. 



E. Comparison with experiment 

In this section, we compare the results of our theoret- 
ical calculations with our experimental results. 

We first compute the near band gap absorption spec- 
tra for a; = and 0.12. This is important since in our 
theoretical calculations, we have only included the effect 
oi sp — d interaction of the Mn ions with the electrons 
and holes. An additional, non-magnetic effect of the Mn 
ions could be to change the band gap with doping x sim- 
ilar to the band gap shift of AlxGai-xAs with increas- 
ing Al content. Figures IT^ a) and ll2r b) show calculated 
near-bandgap absorption spectra for the x ~ Q and 0.12 
samples at (a) 30 K and (b) 290 K without taking into 
account a change in band gap with Mn doping. The 
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FIG. 12: Theoretical absorption coefficient for linear polar- 
ization as a function of photon energy at 30 K (a) and 290 K 
(b) at -B = for Mn concentrations x — 0% (solid lines) and 
X = 12% (dashed lines). Shifts in the bandgaps are entirely 
due to carrier filling effects. These figures should be compared 
with the experimental data shown in Fig. |21 



theory successfully reproduces the experimental results 
shown in Fig. [3 The blue shift of the band gap in the 
a:: = sample is entirely due to band filling effects, i.e., 
the Burstein-Moss shift. The x — sample has a rela- 
tively high electron density. As Mn doping is increased, 
the electron density decreases since the Mn ion acts as an 
acceptor. The theory curves in the figure also show the 
sharpening of the band edge at low temperatures, con- 
sistent with the experiment. The fact that our theory 
curves are able to reproduce the x dependent absorption 
spectra without changing the band gap indicates that any 
X dependence to the band gap is small and we therefore 
neglect it. 

Figures ICT a) and ICT b) show the calculated CR ab- 
sorption coefficient for electron-active circularly polar- 
ized 10.6 /im light in the Faraday configuration as a func- 
tion of magnetic field at 30 K and 290 K, respectively. 
Densities for each sample are given in Table I. In the 
calculation, the curves were broadened based on the mo- 
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FIG. 13: Calculated CR absorption as a function of magnetic 
field at 30 (a) and 290 K (b). The curves were calculated 
based on the Pidgeon-Brown model and the golden rule for 
absorption. The curves were broadened based on the mobili- 
ties reported in Table I. 



bilities of the samples. The broadening used for T = 30 
K was 4 meV for 0%, 40 meV for 2.5%, 40 meV for 5% 
and 80 meV for 12%. For T = 290 K, the broadening 
used was 4 meV for 0%, 80 mcV for 2.5%, 80 meV for 
5% and 80 meV for 12%. 

At r = 30 K, we see a shift in the CR peak as a func- 
tion of doping in agreement with Fig. ^ (a) . For T = 
290 K, we see the presence of two peaks in the pure InAs 
sample. The second peak originates from the thermal 
population of the lowest spin-down Landau level. The 
peak does not shift as much with doping as it did at 
low temperature. This results from the temperature de- 
pendence of the average Mn spin. We believe that the 
Brillouin function used for calculating the average Mn 
spin becomes inadequate at large x and/or high temper- 
ature due to its neglect of Mn-Mn interactions such as 
pairing and clustering. 



V. SUMMARY AND CONCLUSIONS 

We presented a theory for the electronic and magneto- 
optical properties in n-type narrow gap Ini_j;Mn2;As 
magnetic semiconductor alloys in an ultrahigh external 
magnetic field, J5, oriented along [001]. We find several 
key results : i) There is a shift in the cyclotron resonance 
with Mn doping which is not predicted in simple mod- 
els. To lowest order, this shift depends upon (a — f3) and 
can be thought of as a renormalization of the energy gap 
in a magnetic field. The value of the energy gap infiu- 
ences the amount of conduction- valence band mixing, ii) 
Even with no Mn doping, there is spin-splitting in the cy- 
clotron resonance which results from the nonparabolicity 
of the conduction band in the narrow gap material, iii) 
The relative heights of the spin-split cyclotron absorption 
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peaks can allow one to extract information about carrier 
density, iv) At high temperartures and high (> 10%) Mn 
concentrations, the calculated shift of the cyclotron reso- 
nance is not as large as the experimentally observed shift. 
This probably results from the inadequacy of the Bril- 
louin function used for calculating the average Mn spin 
at large x and/or high temperature due to its neglect of 
Mn-Mn interactions such as pairing and clustering 

In modeling the cyclotron resonance experiments we 
used an 8 band Pidgeon-Brown model generalized to in- 
clude the wavevector dependence of the electronic states 
along kz as well as the s-d and p-d exchange inter- 
actions with the localized Mn (i-electrons. Calculated 
conduction-band Landau levels exhibit effective masses 
and g factors that are strongly dependent on tempera- 
ture, magnetic field, Mn concentration x, and k^. At 
low temperatures and high x, the sign of the g factor 
is positive and its magnitude exceeds 100. CR spectra 
are computed using Fermi's golden rule and compared 
with ultrahigh-magnetic-field (> 50 T) CR experiments, 
which show that the electron CR peak shifts sensitively 
with X. Detailed comparison between theory and exper- 
iment allowed us to extract the s-d and p-d exchange 
parameters, a and /3. We showed that not only a but 
also P affects the electron mass because of the strong in- 
terband coupling in this narrow gap semiconductor. In 



addition, we derived analytical expressions for effective 
masses and g factors using the Kane model, which in- 
dicates that {a — (3) is the crucial parameter that de- 
termines the exchange interaction correction to the cy- 
clotron masses. These findings should be useful for de- 
signing novel devices based on ferromagnetic semiconduc- 
tors. 



APPENDIX A: LANDAU HAMILTONIAN 

In this appendix, we write down the Landau contribu- 
tion to the Hamiltonian in the matrix eigenvalue problem 
(|28|l for an arbitrary Landau quantum number n. The 
Landau Hamiltonian matrix can be obtained by taking 
the operator form of the Landau Hamiltonian in Eq. Q 
and operating on the envelope function wavefunction H27() 
making use of the properties of the creation and destruc- 
tion operators and a. The resulting matrix is 



mo 



LV Li 



(Al) 



where the submatrices are given by 



L' = 



-iV'V^^T^ (72 



7i2(2n-3) 



'l23\/3n{n - 1) 



-n23\^Qn{n - 1) 



iV 



-l2Z\J'in{n - 1) 



n 



i723y/Qn{n - 1) 
-(72 + f )fc' - 7i2(2n + 1) iV2{-j2k^ + j'^in + i)) 
-iy2(-72fc2 + 7^'(" + I)) -A' - 71^ - y^in + i) 



(A2) 
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-iv 



n 



-V'^fr 



(72 - f 



7i2(2n + 3) 
l2S.^/in{n + 1) 



-n23\/6n(n -h 1) 



-l23y/3n{n + 1) 



-(72 
-iV2{-j2k' 



«723\/6n(n + 1) 



-7i2(2n-l) zV2(-72fc2+7^'(n-i)) 



+ 7^'(n-i)) -A' 
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L' 






-V 



-«73\/6(^ 



73 



V 



iV 



ilz^/^n- 1) 


-87^ 



-373 V" 





(A4) 



where fc = fc^ is the wavevector along the magnetic field manifold of states. The submatrix L]! is obtained by 
direction, and n is the Landau quantum number for the taking the Hermitian adjoint of L'^. In Eqs. (|A2|) . IjASp . 
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and ljA4p we make the following definitions: 



V 



A 



Here Eg and A are the band gap and spin-orbit splitting 
energies, A is the magnetic length defined in Eq. ((HJ, is 

( A5) the Kane momentum matrix element defined in Eq. Hll|l , 
and the Y^s are the effective mass parameters defined in 
Eq. H12() . The total Hamiltonian, i7„, to be diagonalized 
in the eigenvalue equation H28(l is the sum of the Landau 

(A6) Hamiltonian matrix (|A1(I and the Zeeman and exchange 
Hamiltonians ((TH|l and if^ . 

We note from Eq. ljA4p that the submatrix L'^ is pro- 
(A7) portional to k and so vanishes at A: = 0. In this limit, 
Hn, is block diagonal with respect to the upper and lower 
Bloch basis sets defined in Eqs. Q and 



7,' = f (* = 1...4) 



7i = ^ = 1 . . . 4) 



(AS) 
(A9) 



APPENDIX B: OPTICAL MATRIX ELEMENTS 



(AlO) In this appendix we write down the momentum matrix 

elements used in the computation of optical matrix ele- 
ments in Eq. (|50l) . For the Bloch basis states defined in 
Eqs and (jSJ, the matrix elements for the momentum 

(All) operators Px, Py, and Pz are given by 
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-iV 
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where V is the Kane matrix element defined in Eq. Hll|l . 



ACKNOWLEDGMENTS 






-iV 
-V 







r 



V 











iV 



(B3) 



tion through grant DMR 9817828, and the NEDO Inter- 
national Joint Research Program. 



This work was supported by DARPA through grant 
No. MDA972-00-1-0034, the National Science Founda- 



* Present address: Physics Department, University of Utah, 
Salt Lake City, Utah 84112, U.S.A. 

t Present address; Department of Physics, Faculty of Sci- 
ence, Okayama University, Okayama, Japan. 

^ H. Munekata, H. Ohno, S. von Molnar, A. Segmiiller, L. L. 
Chang, and L. Esaki, Phys. Rev. Lett. 63, 1849 (1989). 

^ H. Ohno, A. Shen, F. Matsukara, A. Oiwa, A. Endo, 
S. Katsumoto, and Y. lye, Appl. Phys. Lett. 69, 363 
(1996). 

^ S. Koshihara, A. Oiwa, M. Hirasawa, S. Katsumoto, Y. lye, 
C. Urano, H. Takagi, and H. Munekata, Phys. Rev. Lett. 
78, 4617 (1997). 

* H. Ohno, D. Chiba, F. Matsukara, T. Omiya, E. Abe, 
T. Dietl, Y. Ohno, and K. Ohtani, Nature 408, 944 (2000). 

^ H. Ohno, H. Munekata, T. Penney, S. von Molnar, and 
L. L. Chang, Phys. Rev. Lett. 68, 2664 (1992). 

® H. Munekata, A. Zaslavski, P. Fumagalli, and R. J. Gam- 
bino, Appl. Phys. Lett. 63, 2929 (1993). 
J. K. Furdyna, J. Appl. Phys. 64, R29 (1988). 

* Y. H. Matsuda, T. Ikaida, M. Miura, S. Kuroda, F. Takano, 
and K. Takita, Phys. Rev. B 65, 115202 (2002). 

^ M. A. Zudov, J. Kono, Y. H. Matsuda, T. Ikaida, N. Miura, 

H. Munekata, G. D. Sanders, Y. Sun, and C. J. Stanton, 
Phys. Rev. B 66, 161307(R) (2002). 

^° K. Nakao, F. Herlach, T. Goto, S. Takeyama, T. Sakak- 
ibara, and N. Miura, J. Phys. E: Sci. lustrum. 18, 1018 
(1985). 

" C. K. Pidgeon and R. N. Brown, Phys. Rev. 146, 575 
(1966). 

J. Kossut, Semiconductors and Semimetals 25, 183 (1988). 
A. L. Efros and M. Rosen, Phys. Rev. B 58, 7120 (1988). 

I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, J. 
Appl. Phys. 89, 5815 (2001). 

J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 



C. Dresselhaus, A. F. Kip, and C. Kittel, Phys. Rev. 98, 
368 (1955). 

" G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

H. Munekata, H. Ohno, S. von Molnar, A. Harwitt, 
A. Segmiiller, and L. L. Chang, J. Vac. Sci. Technol. B 
8, 176 (1990). 

H. Ohno, H. Munekata, S. von Molnar, and L. L. Chang, 
J. Apps. Phys. 69, 6103 (1991). 

S. von Molnar, H. Munekata, H. Ohno, and L. L. Chang, 
J. Magn. Magn. Mater. 93, 356 (1991). 

21 E. O. Kane, J. Phys. Chem. Solids 1, 249 (1957). 

22 J. G. Mavroides, in Optical Properties of Solids, edited by 

F. Abeles (North-Holland Publishing Company, Amster- 
dam, 1972), pp. 351-528. 

2^ W. Zawadzki, in Landau Level Spectroscopy, edited by 

G. Landwehr and E. I. Rashba (North-Holland, Amster- 
dam, 1991), vol. 1, pp. 483-512. 

2* M. A. Kinch and D. D. Buss, J. Phys. Chem. Solids 32, 
461 (1971). 

2^ B. D. McCombe and R. J. Wagner, in Proc. 11th Int. Conf. 
Phys. Semicond. (PWN-Polish Sci. Publ., Warsaw, 1972), 
p. 321. 

2® B. D. McCombe and R. J. Wagner, in Advances in Elec- 
tronics and Electron Physics, edited by L. Marton (Aca- 
demic Press, New Yord, 1975), pp. 1-78. 

2^ L. M. Roth, B. Lax, and S. Swerdling, Phys. Rev. 114, 90 
(1959). 

2* F. Bassani and G. P. Parravicini, Electronic States and 
Optical Transitions in Solids (Pergamon, New York, 1975). 

28 Y. P. Varshni, Physica 34, 149 (1967). 

^° E. D. Palik and R. F. Willis, Phys. Rev. 123, 131 (1961). 
J. I. Pankove, Optical Processes in Semiconductors (Dover, 
New York, 1971). 

^2 C. Rigaux, Semiconductors and Semimetals 25, 229 



20 



(1988). 



